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Abstract. The neutron-proton mass difference in finite nuclei is studied in the framework of a medium- 
modified Skyrme model. The possible interplay between the effective nucleon mass in finite nuclei and the 
Nolen-Schiffer anomaly is discussed. In particular, we find that a correct description of the properties of 
mirror nuclei leads to a stringent restriction of possible modifications of the nucleon's effective mass in 
nuclei. 
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1 Introduction 

An essential quantity in nuclear physics, the effective neu- 
tron-proton mass difference in nuclear matter, Am*, is 
still not known empirically [Tpl On the other hand, there 
exist very different theoretical predictions of this quan- 
tity for isospin-asymmetric nuclear matter [2,3,4,5, GUTUSl 
l??UTnVTTVT^nT?yT^UTK] where the results differ both quali- 
tatively and quantitatively. Such studies of the effective 
neutron-proton mass difference inside nuclei may be rele- 
vant to resolve the Nolen-Schiffer anomaly (NSA) in nu- 
clear physics [TBirrT] and for applications in nuclear astro- 
physics [18, 19, 2 OJ . Although there are many theoretical 
approaches devoted to the explanation of the NSA [21.32, 
|22P5P31|^ll21glgg|[ggi|5n^ this phenomenon is still 
not fully understood. 

In this context, we have recently proposed an effec- 
tive Lagrangian [33] which can be applied to investigate 
isospin-breaking effects in the baryonic sector of the Skyr- 
me model under medium modifications. In ref. [33] the sin- 
gle nucleon properties and Am* lp are studied in infinite, 
asymmetric nuclear matter with a spatially constant den- 
sity, where the surrounding medium influences were taken 
into account as external parameters. Moreover, the model 
can be extended to the evaluation of nucleon properties in 
finite nuclei. Such kind of studies, however with only par- 
tial isospin-splitting effects, have already been performed 
in refs. [34jl35j. It will be interesting to extend these stud- 



1 Thereafter an asterix in an expression indicates that the 
corresponding quantity is explicitly medium-modified. 



ies to the strong and electromagnetic isospin-breaking ef- 
fects in finite nuclei. Consequently, in the present work we 
consider the neutron-proton mass difference in finite nu- 
clei by studying the single-nucleon effective properties and 
by trying to explain the Nolen-Schiffer anomaly within a 
medium-modified Skyrme model. 

The paper is organized as follows: In sect. [5] we re- 
view the model Lagrangian which was adapted in ref. |33j 
to skyrmions in asymmetric nuclear matter. Furthermore, 
we specify the input parameters and the medium function- 
als, which express the influence of the surrounding envi- 
ronment onto the single-nucleon properties, together with 
the parameterization of the neutron- and proton-density 
distributions. Section [3] centers on the spherically sym- 
metric case where a classical skyrmion is located at the 
center of a spherical nucleus. In particular, we contrast 
the exact solution of the equations of motion with a vari- 
ational approximation and compare in sect. 3] the results 
of both calculations. Next, in sect. [5] we extend our stud- 
ies to the nonspherical case where the nucleons are lo- 
cated at some finite distance from the center of nucleus. 
We discuss the consequences of the nonspherical case for 
the time-dependent Lagrangian, the appropriate classical 
solutions and the quantization procedure. In sect. we 
present the pertinent calculations for the electromagnetic 
part of the neutron-proton mass difference in finite nuclei 
and in sect. [7] we discuss the results of the nonspheri- 
cal approximation. Our discussions of the Nolen-Schiffer 
anomaly are presented in sect. [5] and a summary of the re- 
sults is given in sect. [HI Finally and for completeness, two 
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appendixes are added which elaborate on the nonspherical 
scenario: the first is devoted to the mass functional and 
moments of inertia (appendix |X| and the second to the 
pertinent charges ( appendix [B| . 



2 Lagrangian of the model 

We start with the Skyrmc-modcl Lagrangian presented in 
ref. [33] 
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where Einstein's summation convention is always assumed 
(if not specified otherwise). The chiral SU{2) matrix U 
has the form U = exp(2ir a 7r a /i 7 ' 7r ), where 7r a (a = 1, 2, 3) 
are the Cartesian isospin-components of the pion field. 
F w = 2/, is the pion-decay constant, while e is the di- 

mensionless Skyrme constant. A4- = ^/(m^ ± — m 2 r0 )/2 

is defined in terms of the masses of the charged and neu- 
tral pions, where for the sake of convenience the mass 
of the neutral pion is denoted as m, = m^o in ([4]). (|6j) 
and flU). Finally, x°°, xf, Xp, A Xs , and A Xp represent 
the medium functionals which express the influence of the 
surrounding environment onto the singlc-skyrmion prop- 
erties. The numerical values of the input parameters in the 
present work have slightly been changed in comparison to 
the work |33j . The parameters F^ and e are fixed in such 
a way that the free-space (density p = 0) PDG-valucs 
of the proton and neutron masses (m p = 938.27 MeV 
and m n = 939.56 MeV [36]) together with the empirical 
(isospin- averaged) mass of the delta (rriA = 1232 MeV) 
are reproduceq 2 ]. Furthermore, as input for the free mass 
of the neutral pion we still use the PDG- value 36J: m T o = 



Note that in the previous work [33] the input for F, and e 
was determined by a fit to the isospin-averaged masses of both, 
the nucleon and the A. 



134.977 MeV. All of these choices together induce the fol- 
lowing values for the input-parameters: F^ = 108.783 MeV, 
e = 4.854, and m w ± = 135.015 MeV; the latter is the 
strong part of the mass of the charged pion. Note that the 
dominant electromagnetic contribution to m^± — m^o is 
still beyond the scope of the model. 

In the Lagrangian ([I])-© the term £* ym expresses the 
isospin-symmetric part, whereas AC mcs and AC* nv are the 
isospin-breaking terms that arise from the explicit symme- 
try breaking in the mesonic sector and the isospin asym- 
metry of the surrounding environment, respectively. 



2.1 Medium functionals 



The medium functionals are given as in 
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where 



P{r) = p n (r) +p P (r), 
6p(r) = p„(r) -p p (r), 
X] = 1 + m^/m^- 



(9) 
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Here k-p = [3tt 2 p/2} 1 / 3 is the total Fermi momentum, 
Pnir) and p p (r) are the neutron- and proton-distribution 
densities in the nucleus, mn is the isospin-averaged nu- 
cleon mass: ttin = (m p + m n )/2. The parameters of the 
medium functionals have the values bo = — 1.206m 7i r 1 , 
61 = -1.115m" 1 , c = 0.21m" 3 , ci = 0.165m" 3 and 
g' = 0.47 [33] ■ The relations between the parameters 60,1 
and the s-wave pion nucleon scattering lengths 6 ,i, based 
on chiral perturbation theory at order C(m 3 ), can be 
found in refs. [571138] . As input for the p-wave scatter- 
ing volumes, Cq and ci, we use the threshold values of the 
'current' SAID analysis [Hof 5 ! . 



2.2 Proton and neutron distributions in finite nuclei 

From now on we will only concentrate on nuclei in the 
ground state and furthermore consider only nuclei which 
either are "magic" or which are very near to shell closure. 



3 For more references and explanations about the chosen val- 
ues of 60,1 and co,i see ref. [33] . 
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This allows us to use spherically symmetric approxima- 
tions to the nucleon distributions inside the nuclei. Conse- 
quently, the distribution densities of protons and neutrons 
in finite nuclei are parameterized in the standard way [40] 
as 



PpA r) = (Z,A-Z) I - Jr h 



p. 11 
1 



9 fl 
* "p,n 



p,n 



l + exp{(r-r' )/a' } ' 



(18) 



in terms of the mass-number- and isospin-dependent pa- 
rameters 



r (o) A i/3 

.(1) 
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' p,n' v ; 



~p,n 



(19) 



In these expressions A is the mass number, Z is the num- 
ber of protons, and A = (A—2Z) j A is the isospin-asymmet- 
ry parameter. As input for the parameters of the neutron- 
and proton-distributions we take the values 



r (°) 

r (D 
'p 

-p 2) 
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„(2) 



1.2490 fm, 
-0.5401 fm, 
-0.9582 fm, 

0.4899 fm, 
-0.1236 fm, 



.(0) 



.(2) 
n 

(2) 



1.2131 fm, 
-0.4415 fm, 
0.8931 fm, 
0.4686 fm, 
0.0741 fm. 



These values were obtained from calculations based on 
effective Skyrme interactions used in the description of 
the properties of finite nuclei [41] . 

As our main interest is the Nolen-Schiffer anomaly, we 
will concentrate our investigations to the following pairs of 
mirror nuclei: 15 N versus 15 0, 17 versus 17 F, 39 K versus 
39 Ca, and finally 41 Ca versus 41 Sc. An important point in 
the calculation of the effective nucleon properties is the 
proper definition of the nuclear background density. For 
example, by adding cither a neutron or a proton to the 
background of the spherical nucleus 16 O one can get the 
nucleus 17 or 17 F, respectively. We assume that possi- 
ble changes in the structure of the 16 O core due to the 
additional valence nucleon are small and can be ignored. 



3 Parameterization of the classical solutions I 
(spherically symmetric configurations) 

When the skyrmion is located at the center of the nu- 
cleus, one can still use the spherically symmetric hedgehog 
ansatz 

U = exp[ir • (r/r)F(r)} (20) 

and, following the two-step method of refs. [33II42] . con- 
struct the time-dependent Lagrangian in terms of the stan- 
dard angular velocities u>i of the collective modes and the 
constrained angular velocity a* . Then the time-dependent 



Lagrangian is given as 



£*d 3 r = -Mt v -MtA_ + ^-A* 



-oj 3 (a*A* + A* 



■ I — A* + A 



(21) 



where 
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(22) 



+ Fl (ml + X °°) (l-cosF)|r 2 dr, 



is the in-medium mass of the solution when it is not per- 
turbed (NP) by any isospin-breaking effects in the mesonic 
sector or the nuclear environment. The abbreviations F r = 
dF/dr and S = smF have been used, where F = F(r) is 
the chiral profile function of the hedgehog ansatz. 



2tt 

A- = —F. 



S z r z dr 



(23) 



is a moment-of-incrtia-type term resulting from the ex- 
plicit (strong) isospin-breaking due to the pion masses, 
whereas 



oo 

A* = ?pJ j (l + m~\f)S^dr 



(24) 



is the moment of inertia of the skyrmion in the nuclear 
medium. Finally, 



A* = ^F*J AaS 2 r 2 dr 
o 



with 



Aa 



2m 7I 



(A Xs + A Xp ) 



(25) 



(26) 



is the isospin-breaking factor due to the isospin asymme- 
try of the nuclear environment. We note that on top of 
the terms considered in the previous refs. [3^114^] also den- 
sity gradients of the surrounding environment are relevant 
which result from the p-wave pion-nucleus scattering. 

In order to recover the minimization functional in the 
classical approximation (cj, = 0), one applies the con- 
straint [42] 

a* 2 = 2M 2 _A_/A* (27) 
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such that from eq. (|21|) the Lagrangian 

"A* 



L* = —Mi 



NP 



(28) 

is generated for vanishing a;,. The variation of the La- 
grangian (|28|) gives the pertinent equation of motion for 
the hedgehog profile function F(r), 

F%(1 - X ° P ) (r 2 F rr + 2rF r - S 2 ) - F 2 X ° p , r r 2 F r 
4 



(29) 



-F 2 (ml + X °°) Sr 2 + ^UaS 2 r 2 = , 



where the additional abbreviations F rr — d 2 F/dr 2 , S2 = 
sh^i* 1 and Xpr = dXp/dr were introduced. The boundary 
conditions 



lim F(r) = tt — Cr , 
lim F(r) =D (1 + m^r) 



(30) 



where C and D are constants, correspond to a classical 
soliton of baryon number B = 1. 

As an alternative, we will also use the following trial 
function for approximating the exact solutions of the above- 
given differential equation (f29l) : 



F(r) = 2 arctan 1^(1 + m n r) \ e - f{r)r , 



(31) 



Here r$, /3o, Pi, and /?2 are variational parameters. Note 
that our trial function F(r) differs from the widely used 
trial function 2 arctan {r 2 /r 2 }, since it takes into account 
not only the asymptotical behavior at the origin and the 
far-distance region, but approximates also the behavior 
in the intermediate range. The trial ansatz will be useful 
in our following studies of approximate solutions in the 
nonsphcrical case. 

As in the previous work [53] the final value of a* is 
found by an iteration procedure, starting at a* = 0. 

Whenever the skyrmion can be parameterized as a 
spherically symmetric hedgehog configuration, the quanti- 
zation procedure and the expressions for the electromag- 
netic part of the neutron-proton mass difference remain 
the same as given in in ref. |33J, except for the fact that 
some global medium functionals acquire a local density 
dependence (see eq.s ©-(HSJ)) and that the expression of 
A* also has an additional contribution from p-wave pion- 
nuclcus scattering]. 



4 Results in the spherically symmetric 
approximation 

The results arising from (i) the direct solutions of equa- 
tions of motion (|29[) and (ii) from the approximate varia- 
tional procedure with the trial function (f3"Tj) are presented 



4 Compare the present expressions 
expressions (40) and (38) of ref. [33 j - 
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Fig. 1. The solutions of eq. (|29|) . The solid curve represents the 
profile function of a skyrmion in free space, the dashed curve 
stands for the profile function of a skyrmion located at the 
center of an le O background, and the dotted curve represents 
the case of a skyrmion located at the center of a 40 Ca core. The 
approximate trial function (|31[) gives almost identical results 
(indistinguishable at the scale of the figure) for each of these 
three cases. 



in Table [T] for several nuclei. One can see that the trial 
function (|3ip very well approximates the actual solutions 
of the model. The deviations from the exactly calculated 
static properties of the nucleon due to the approxima- 
tion (f3~T|) are in the range of 0.01 -j- 5.4 %. This can be 
seen by comparing the lines (i) and (ii) for each nucleus 
presented in Tabic [TJ In free space the deviations are even 
smaller: the maximum deviation is in Am np and is ap- 
proximately equal to 2%. 

As expected, in accordance to our previous calcula- 
tions for infinite nuclear matter in the isospin symmetric 
case [33], the neutron- proton mass difference is slightly 
increased in comparison to the free space case. Also note 
that in this model the effective nucleon mass has a strongly 
decreased value at the center of nucleus. We will return to 
this point later on. 

For illustrative purposes the solutions of the differen- 
tial equation (|29p for a skyrmion in free space as well as 
inside the nuclei 16 O and 40 Ca are shown in fig. [T] Note 
that the profile function is rather insensitive to small den- 
sity variations: the solutions of the nuclei 14 N and 16 or 
38 K and 40 Ca, respectively, almost coincide. Finally, from 
the results presented in Table [T] it is clear that the trial 
function (|3ip with appropriately chosen parameters is al- 
most as good as the true solution. For this reason the trial 
functions are not presented in fig.HJ since graphically they 
cannot be distinguished from their exact counterparts. 



5 Nonspherical solutions 

5.1 Time— dependent Lagrangian 

When the skyrmion is located at some finite distance from 
the center of the nucleus, the spherically symmetric hedge- 
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Table 1. The values of the variational parameters and the static properties of nucleons when the pertinent skyrmion is either 
in free space or added to the center of a finite (background) nucleus. Here m* is the in-medium proton mass, Am^ p is the 
in- medium neutron-proton mass difference and Zimn P EM ' is its electromagnetic part; /i p and /i,* are the in- medium proton and 
neutron magnetic moments in units of nucleon Bohr magnetons (n.m.) in free space; {r 2 }*^~g 2 and (r 2 )*^ 2 are the in-medium 
isoscalar (S) and isovector (V) charge radii of the nucleons. For a given nucleus we present (i) the results corresponding to 
the exact solution of the differential equation (|29p and (ii) the results of the approximate variational procedure where the trial 
function (|31[) has been used. 



Element 


ro 
[fm] 


10/3 
[ma] 


Pi 
[m w ] 


P2 

[ml] 


m* 

[MeV] 


Am* p 
[MeV] 


A *(BM) 

[MeV] 


[n.m.] 


Mi* 

[n.m.] 


[fm] 


(r%% 2 
[fm] 


free 
space 


(i) 
(ii) 


0.954 


0.075 


1.311 


-0.009 


938.268 
938.809 


1.291 
1.313 


-0.686 
-0.687 


1.963 
1.966 


-1.236 
-1.241 


0.481 
0.481 


0.739 
0.739 


14 N 


(i) 
(ii) 


1.393 


0.076 


0.920 


0.226 


593.285 
598.505 


1.668 
1.655 


-0.526 
-0.536 


2.355 
2.230 


-1.276 
-1.209 


0.656 
0.648 


0.850 
0.810 


16 Q 


(i) 
(ii) 


1.426 


0.076 


0.907 


0.219 


585.487 
590.175 


1.697 
1.685 


-0.517 
-0.527 


2.393 
2.341 


-1.297 
-1.232 


0.667 
0.660 


0.863 
0.825 


38 K 


(i) 
(ii) 


1.493 


0.076 


0.841 


0.153 


558.088 
559.957 


1.804 
1.802 


-0.480 
-0.485 


2.584 
2.550 


-1.422 
-1.377 


0.722 
0.718 


0.942 
0.910 


40 Ca 


(i) 
(") 


1.489 


0.076 


0.839 


0.149 


557.621 
559.378 


1.804 
1.802 


-0.478 
-0.483 


2.569 
2.557 


-1.428 
-1.383 


0.724 
0.720 


0.947 
0.914 



hog ansatz cannot be used anymore, since the background 
profile is not spherically symmetric. Thus there exist fur- 
ther deformations in the isotopic as well as in the coor- 
dinate space j34ll35j . In this case one should use a more 
generalized ansatz for the matrix U : 

U(r-R)=exp[iT-N{r-R)P(r-R)] , (32) 

where R is the distance between the geometrical centers 
of the skyrmion and the nucleus (see fig.[2|). Furthermore 
note that the argument of the profile function P depends 
on the direction in coordinate spac<|3. For consistency, the 
isotopic vector TV must also be generalized as 

(sin6>(? — R) cos^X 
sm0(r - R) sin ip (33) 
cos0(r-R) J 

in terms of a distance- and direction-dependent general- 
ized angle = 0(r — R). Furthermore, due to the az- 
imuthal symmetry of the problem (again sec fig. [2|) , one 
can restrict the form of both functions as 

P = P(\r-R\,6), 9 = 0(\r-R\,9) (34) 

and put the origin of the coordinate system at the geo- 
metrical center of the skyrmion: r' = r — R. For nota- 
tional convenience we rename r' to r, keeping in mind 
that the medium functionals acquire an R dependence 
(see appendix [Ajf^l . 

5 Our notation distinguishes between spherically symmetric 
profile functions (F) and direction-dependent ones (P). 

6 This redefinition procedure is explained in more detail in 
ref. [34]. 




Fig. 2. A sketch of a skyrmion located inside a finite nucleus. 
Here R is the distance between the center of the nucleus (O) 
and the geometrical center of the skyrmion (O'); the vectors 
(angles) r (9) and r' (9') refer to the body-fixed coordinates of 
the nucleus and skyrmion, respectively. Without loss of gener- 
ality both coordinate systems can be orientated in such a way 
that their z-axes coincide. 



In order to apply the collective-coordinatc-method to 
the configuration given in (f32l) - (134l) we consider the fol- 
lowing independent rotations in coordinate space and in 
isotopic space 

P = P(K- 1 (t)r), N =I{t)N (K- 1 ^) , (35) 
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where TZ and X are rotational and iso-rotational matri- 
ces, respectively. Applying these time-dependent rotations 
to the two-step method of ref. [33] , one arrives at a La- 
grangian which - after a spatial integration - is given as 

T* - A f* M 2 A I W 1 +UJ 2 a* 

Q 2 + Q% 

- (wi«i+w 2 i?2)vl* fiil2 + 2 A* an>12 

i 2 7i wi?,33 

+ (u> 3 -n 3 + a*)At nv . (36) 

Here uji and i?; arc the angular velocities in isotopic and 
coordinate space, respectively. The explicit expressions for 
the functionals, M^ P [P, 0] and each A[P, 0], can be found 
in appendix IA1 



where the {7™} are variational parameters and the cosine 
functions are chosen to maintain periodicity in 9. Simi- 
larly, £ can be chosen as 

oo 

C(r, 6>) = re- s « r2 ^ $n sin 2n6 , (41) 

71=1 

where the {S n } are also variational parameters. Note that 
the arguments of the sine-functions are chosen to be a 
multiple of 29 in order to avoid singularities of the form 
sin 0/ sin 9 in the Lagrange functional (see appendixlAl 
for the explicit form of the Lagrange functional) . Further- 
more, the r dependence of £ is chosen in such a way that 
the equalities 0(0,9) = 9 and 0(oo,9) = 9 are repro- 
duced. 

Finally, in terms of the trial functions |39|) , (|40)l and 
(|4Tj) . the Lagrange functional (f38|) will be extremized. 



5.2 Parameterization of the classical solutions II 



5.3 Quantization of the nonspherical solutions 



As in our previous works [33U42] . by incorporating the 
constraint for the angular- velocity parameter a* , which in 
the azimuthally symmetric case has the form 



a* 2 = 2M 2 



A, 



A 



(37) 



w Q, 33 



one arrives at the following Lagrange functional at the 
classical level: 



L* = -M* P + a* At 



(38) 



In order to extremize the Lagrange functional (|38[) . one 
would have to solve the complicated system of coupled 
partial differential equations that arises from the variation 
of the Lagrangian. To avoid all the unnecessary technical 
difficulties connected with a numerical solution of that 
equation system, we rather use an approximate variational 
procedure that is specified as follows: first of all, as fig. Q] 
indicates, the profile function is only weakly dependent on 
the density of the medium (compare the solid and dashed 
curves of fig. [J) and almost insensitive to small density 
changes (compare the dashed and dotted curves of fig- LXI) - 
Consequently, we can apply the approximations 



P(r,9) = 2arctan<^ +m^r)(l + u{6)) \ e - f ^ r 



0(r,6) = 9 + ((r,6). 



(39) 



where the special case P(r, 0) = F(r) is the spherically 
symmetric profile function (in the form as presented in 
eq. (|3ip ) and where the functions u and £ satisfy the in- 
equalities |u| < 1 and |e| < 1 in the regions rC [0, oo) and 
9 G [0, 71"]. Following the ideas of ref. [33] we use for the 
function u the following ansatz 



Defining canonical conjugate variables in the body-fixed 
reference system 



dL* dL* 
U = -= — and J, = — — 

OUJj Oil; 



(42) 



one obtains from the time-dependent Lagrangian (|36[) the 
Hamiltonian: 



H = M£ P + M 2 _A mes 



A* 



'NP 

^2 i rfi2\ A* 



2A 



tuf2, 33 
2 i r2\ a* 



^(AZu^Aqo^u ~ A*J n i 2 ) 

A* A* — A* 2 

( A* \ - 



T 2 



2.1 



uQ ,33 



.1 



(43) 



wf2,33 



By sandwiching the Hamiltonian between the appropri- 
ate baryon states |T, T3; J, J3 = — T3) one determines the 
energy of a nucleoli inside a nucleus as 



e = m: 



M 2 _A mes + 



A 



*2 
env 



2.1 



to f 2, 33 



^T + fr A2 2 5 n '^ (T(T + 1) - T 2 ) 



T 2 
1 3 



A* 



2.1 



U)Q.33 



A* 



(44) 



i?,33 , 



Consequently, the strong part of the neutron-proton mass 
difference in the interior of a nucleus takes the form 



u(9) 



^] In COS T 

71=1 



(40) 



zim *(strong) 



A* 



A 



(45) 



to f 2, 33 
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6 Electromagnetic part of the neutron-proton 
mass difference in nonspherically deformed 
states 

The electric (E) and magnetic (M) form factors of the 
nucleon are defined through the expressions 

Ge(Q 2 ) = 1 /d 3 re*V(r), 

G* M {q 2 ) = ^ J d 3 re"Hrx i(r)] , (46) 

where q 2 is the squared momentum transfer. Furthermore, 
j° and j correspond to the time and space components of 
the properly normalized sum of the baryonic current B* 
and the third component of the isovector current V* of 
the Skyrme model. 

For the problem at hand, it is advantageous to expand 
the plane wave factor in ggClas 

OO 

e iqr = e iqrcose =J2 il ( 21 + l)^(cos0)Ji(gr) , (47) 

;=o 

in terms of Legendre polynomials Pi and spherical Bessel 
functions <//. In this way we get the final expressions for 
the electromagnetic form factors 

OO 

g:V) = E^ 2Z + 1 ) g ^ 2 ). (48) 

Here 

G%itf) = Jd\p b a (r,e)P l (co S e)J l (qr) (49) 

are partial form factors, where the label a stands either for 
isoscalar (S) or isovector (V) form factors, and the label b 
denotes either electric (E) or magnetic (M) form factors. 
Explicit expressions of the corresponding charge (pg) and 
magnetic (pm) densities are given in appendix [Bj As ac- 
tual calculations show, the absolute values of the partial 
form factors decrease quickly with increasing partial wave 
number I within the present approach. 
Finally, applying the formula |43j 

OO 

Am:^ = -^Jd q {Gl?( q *)GV* {q *) 
o 

-^ G ^ {q2)G ^ ( - q2) } (50) 

where a « 1/137 is the electromagnetic fine structure 
constant, one can calculate the medium-dependent elec- 
tromagnetic part of the neutron-proton mass difference. 

7 We can rotate the total system, background nucleus and 
skyrmion, always such that the z-axis of the reference frame 
coincides with the z-direction of the coordinate system of the 
skyrmion in the body-fixed frame. 



7 Results in the nonspherical approximation 

First of all let us remark that numerically the transitions 
lim <9(r, 0) -> and lim P(r, 0) -> F(r) 

R->0 R->0 

are smooth when the nucleon 'moves' to the center of the 
nucleus, and that in this case the exact solution is repro- 
duced with high accuracy, see Table [1] One might there- 
fore expect that the trial functions (|39p - (|4"Tj) are not too 
far from the the true solutions even when the distance 
between the center of the skyrmion and the center of the 
nucleus is not equal to zero, R =/= 0. 

Obviously, the variational parameters ro, {/?«}, {in} 
and {<5„} depend on the distance R between the centers 
of the skyrmion and the nucleus. Note that the strength 
of the parameter 7„ decreases with increasing index n. 
For example, 73 is one order magnitude smaller than 71 
at those distances where the deviations of the skyrmion 
from the spherical form is maximal. The strength of 8 n de- 
creases even more rapidly with increasing n: 82 is one or- 
der of magnitude smaller than 8\ at the above-mentioned 
distances and £3 is almost zero for all values of R. For 
performing high accuracy calculations it is therefore suffi- 
cient to retain the 7„ parameters up-to-and-including 74 
and the 8 n parameters up-to-and-including 82- 

The behavior of the static properties of the nucleons in 
nuclei and the corresponding variational parameters are 
discussed in a more detail in ref. |34j . Qualitatively, we 
have similar results in the present approach. For example, 
the in-medium mass of the proton starts with the effective 
value presented in Tablc[T]whcn it is near the center of the 
nucleus. Then it increases monotonically with increasing 
distance between its center and the center of the nucleus 
until it approaches the free space value at the border of 
the nucleus, see fig. [3] 

We will come back to the discussion of the effective nu- 
cleon mass in sect. [8] Here, however, we concentrate on the 
effective in-medium neutron-proton mass difference. The 
behavior of Z\mnp strong ^ inside several nuclei is presented 
in fig. [4j The strong part of the effective mass difference 
has a non-monotonic behavior. This is a consequence of 
the density now being a local quantity and of the appear- 
ance of additional isospin-breaking contributions due to 
the density gradients arising from the p-wave pion-nuclcus 
scattering (see eq. ((26]) h One can explicitly see in fig. [4] 
that at the surface of the nucleus, where the density gra- 
dients are large and the local isospin asymmetry is high, 
the value of Z\mn p stIons ^ is at an extremum. 

In fig.[5]thc electromagnetic part of the neutron-proton 
mass difference is shown for several nuclei. The variations 
in the electromagnetic part of the neutron-proton mass 
differences are small compared with their strong counter- 
parts. Qualitatively, their behavior is the same for all nu- 
clei as in the case of strong part. Quantitatively, both 
parts, Zlmnp' 1 ™ 8 ' and An%np , are mainly governed by 
the (local) behavior of the total nuclear density and there- 
fore differ if the total densities of the investigated nuclei 
are different. 
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R[fm] 



Fig. 3. The dependence of the effective mass of the proton 
rrip on the distance 7? between the centers of the skyrmion 
and the nucleus. The solid curve represents the result for an 
additional skyrmion embedded into a 14 N core, the dot-dashed 
curve stands for the case of 16 O, the dashed curve refers to the 
case of 38 K and dotted curve refers to the case of 40 Ca. The 
horizontal line indicates the pertinent free space value of the 
proton mass, see line (ii) in Table [1] i.e. m p = 938.81 MeV. 



2.4 t 



> 




2 4 6 

R[fm] 



Fig. 4. The dependence of Arrinp Ton on the distance R be- 
tween the centers of the skyrmion and the nucleus. The nota- 
tions are same as in fig. [3] The horizontal line indicates the 
free space value of Zimi S p rons ' = 2 MeV (see the corresponding 
line (ii) in Tabled}. 

Finally, in fig. [6] the dependence of the total neutron- 
proton mass difference Am^ p on the distance R is shown 
for several nuclei. 



8 Calculation of the Nolen-Schiffer anomaly 

A long standing problem in nuclear physics is the Nolen- 
Schiffcr anomaly (NSA) observed in mirror nuclei p~6|ri7] . 
Let us consider this point in more detail and show how 
the NSA can be treated in the present framework. 



-0.4 t 




2 4 6 8 



R [fm] 

Fig. 5. The dependence of zimn p EM ' on the distance R between 
the centers of the skyrmion and the nucleus. The notations are 
same as in fig. [3] The horizontal line indicates the free space 
value of Zimip M ' = —0.69 MeV (see the corresponding line (ii) 
in Table [TJ. 




R[fm] 



Fig. 6. The total neutron-proton mass difference Am^ p in 
nuclei. The notations are same as in fig. [3] The horizontal line 
indicates the free space value of Am np — 1.31 MeV (see the 
corresponding line (ii) in Table [T}. 

The mass difference between mirror nuclei which differ 
by one unit in their charges, Z\ — Z<i = 1, is usually split 
into two contributions 

AM = i +1 M N - iM N+1 = AE EM - Am* np , (51) 

where AEem is the computed Coulomb energy difference, 
including various corrections (by, e.g., the exchange term, 
the center-of-mass motion, finite size effects of the proton 
and neutron charges, magnetic interactions, vacuum ef- 
fects, the dynamical effect of the neutron-proton mass dif- 
ference, and short-range two-body correlations). Further- 
more, the quantity Am^ p is the effective neutron-proton 
mass difference. The measurement of the mass difference 
AM and the theoretical calculation of AE-em can be done 
with great accuracy (within 1% error) [17] . Furthermore 
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if one assumes that ^to* p is constant and equal to its free 
space value, then eq. (|51j) will not be satisfied. This phe- 
nomenon is called the Nolcn-Schiffcr anomaly, persistent 
throughout the periodic table. Quantitatively, the NSA 
ranges between few hundred keV for the lightest nuclei 
up to one MeV or more for the heaviest nuclei. A possible 
resolution of the anomaly is the assumption that the effec- 
tive neutron-proton mass difference ^to* p would decrease 
with increasing mass number A. Consequently, the NSA 
would be z^nsa = Am ap — Zim* p . 

Within the present approach Z\m* p has a local R de- 
pendence according to the location of the nucleons inside 
the nuclei as represented in fig. [6] In order to compare 
our results with the experimental data one has to average 
the value of Am^ p with respect to the separation R. In 
this respect we note that the nucleons inducing the Nolen- 
Schiffer anomaly are valence nucleons which, if the mir- 
ror nuclei differ by one particle or hole from the (magic) 
closed-shell case, must be located in the peripheral region 
of each of the mirror nuclei. For example, in a nuclear 
shell model with a Wood-Saxon potential the wave func- 
tions of the valence nucleon of a nucleus of mass number 
A = A mag i C ± 1, respectively, can be fitted very well by 
the Gaussian form 



V>i n ' p) 



(R) 



27rr(3/2 + n A ) 



x cxp ■ 



Rf 



(52) 



where r is the Gamma function and the parameters b A 
and tia have the following valued 



b (p) 


= 0.87177V, 


An) 
°15 


= 0.886m OT , 


rii 5 


= 1, 


b (p) 
°17 


= 0.843777^, 


An) 
"17 


= 0.860777,,, 


nn 


= 2, 


b (p) 
°39 


= 0.785777^, 


°39 


= 0.804m T , 


n 3 9 


= 2, 


l(p) 
°41 


= 0.774777^, 




= 0.794m w , 


7741 


= 3 



The averaged neutron-proton-mass difference is given 



Z\m* p = 777* - 777* , (53) 



m* n jR)(^(R)) 2 d 3 R. (54) 



Defining the difference of neutron-proton probability den- 
sities as 

(^(R)f - {^(R)f ^ A^(R) 
one can rewrite eq. (|53p as 

Am* np * J (Aii 2 J m ;+(^) 2 Am: p ) d 3 i? 

ee Am^ + Am*™ , (55) 

8 For the values of the Wood-Saxon potential parameters see 
ref. [44]. 

9 Note that the wave functions (|52[) are properly normalized. 



ignoring the subleading contribution of the cross term 
jA^Am* np d 3 R. 

Then the Nolcn-Schiffer anomaly for a given pair of 
mirror nuclei with the mass number A simply takes the 
form 



A 



NSA 



ArrfA? + ArrfW 



(56) 



The dependence of the Nolcn-Schiffer anomaly on the 
mass number A is presented in Table [2] (see the third col- 
umn 'a rcn =0, present approach'). One can see that our 
predictions of the behavior of the NSA qualitatively go 
into the right direction. But quantitatively, the results 
are too big by one order of magnitude. This is mainly 
due to the strong shift of Amp (see the third column 
of Table \2§ that results for three reasons: (i) the rather 
large renormalization of the effective nucleon mass, (ii) 
the pronounced R dependence of to* inside the nucleus 
(see fig. and (iii) the relative swelling of the proton 
distributions due to the Coulomb factor, i.e. Aip np ' ^ 0. 
For example, the averaged in-medium mass of the valence 
proton in 17 is reduced to 777* = 812.35 MeV. This drop 
of about 125 MeV is very large in comparison with the 
empirical value of the binding energy per nucleon in nu- 
clear matter. For heavier nuclei, where the density in the 
interior approximates the normal nuclear matter density 
(sec fig. [3]), the drop of the averaged effective mass is even 
larger, e.g. m p — to* ~ 150 - 200 MeV around 40 Ca (see 
the second column of Table[2]) down to ~ 300 MeV around 

208p b _ 

If only the contribution Amnp 2 ^ (due to the explicit R 
dependence of the neutron-proton mass difference) were 
considered, then the NSA in the present approach would 



np 



ArfC^ 



< 0. 



even have a negative sign: Am 

In our model the effective mass of the nucleon strongly 
depends on the phcnomcnological input parameters of the 
pion- nucleus sector, mainly via the p-wave scattering vol- 
ume Co and the Landau parameter g' (see subsection l2.1[) . 
In ref. [45] , according to the systematics of the effective ax- 
ial coupling g* A in finite nuclei of the range 5 < A < 39 46 1 , 
the values cq = 0.15to~ 3 and g' = 0.6 were predicted^. 
These values imply a smaller drop of the effective mass 
of the nucleon, e.g. to*^/ton ~ 0.76 in normal nuclear 
matter [45]. Although, the use of smaller values of c 
and bigger values of g' improves the NSA results {e.g., 
Z NSA = 2.86 MeV for the pair 15 0- 15 N for c = 0.15m^ 3 
and g' = 0.75), this is still not enough to satisfy the ex- 
perimental data. 

Instead of fine-tuning these parameters to unphysical 
values, one might rather invert the problem and try to 
estimate the effective nucleon mass inside finite nuclei ac- 
cording to the results in the isospin-breaking sector. To 
perform this task we fine-tune an artificially added renor- 
malization parameter ct ren in the expression 

TO* p (i?,a rC n) = 777* p (R) + (m n!P - to* p (i?))a ron (57) 



10 Compare with the values Co = 0.21m,,. 3 and g' — 0.47 used 
in the present work. 
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Table 2. The averaged mass of the valence proton in a given nucleus m*, the contributions to the effective neutron-proton 
mass difference (see eq. $55} ) and the comparison of the Nolen-Schiffer discrepancy Z\nsa (see eq.s (156[) , JSTJl) calculated in the 
present approach with other "empirical" results. All quantities are in units of MeV. 



Present approach 



Nuclei 


m* 


Qf rcn = 


a rcn = 0.95 


^NSA 


^NSA 




Qrcn = 


Qrcn = 0.95 


Am^ 


Am*™ 


Ansa 


Am™ 


Ami™ 


ZiNSA 


ref. QI] 


ref. [Tr] 


15 0- 15 N 


767.45 


928.30 


-4.27 


1.56 


4.02 


-0.21 


1.33 


0.20 




0.16 ±0.04 


17 F- 17 


812.35 


930.54 


-5.53 


1.52 


5.33 


-0.28 


1.32 


0.27 


0.31 


0.31 ±0.04 


39 Ca- 39 K 


724.78 


926.16 


-8.11 


1.67 


7.75 


-0.41 


1.33 


0.37 




0.22 ±0.08 


41 Sc- 41 Ca 


771.71 


928.51 


-9.74 


1.62 


9.44 


-0.49 


1.33 


0.47 


0.62 


0.59 ±0.08 



of the effective nucleoli mass in such a way that the NSA 
is satisfied. The results are presented in Table [2] (see the 
column 'a ron = 0.95'). It can be seen that a success- 
ful description of the correct order of the NSA implies 
a rather small drop of the mass of the valence nucleons: 
m n>p — m* (a rC n = 0.95) ~ 10 MeV which is close to the 
empirical binding energy per nucleon in nuclear matter. 
In this case the contribution to the NSA from Arn^p 2 ^ can 
be neglected: Am np ~ Am^p (a rcn = 0.95) 0.02 MeV. 



9 Conclusions and Outlook 

In summary, we have studied the effective neutron-proton 
mass difference Am^ p in finite nuclei in the framework of 
a medium-modified Skyrme model. The value of Am^ p for 
a given nucleus approaches an extremum near the nuclear 
surface where the local density gradients are large and the 
isospin asymmetry is high. 

We have discussed the relevance of our results for the 
Nolen-Schiffer anomaly. Qualitatively, our approach pre- 
dicts the correct behavior of the Nolen-Schiffer anomaly. 
But quantitatively it is not satisfactory: the results are one 
order of magnitude too large. Clearly, the part of our cal- 
culation relevant to the Nolen-Schiffer anomaly depends 
on the proton and neutron distributions of the mirror nu- 
clei and is very sensitive to the behavior of the wave func- 
tions of the valence nucleons in the peripheral region of 
the nucleus. We have pointed out the possibility that the 
Nolen-Schiffer anomaly may rather follow from the be- 
havior of the effective nucleon mass in finite nuclei than 
from the effective neutron-proton mass difference: our cal- 
culations imply that the Nolen-Schiffer anomaly could not 

and, maybe, should not be saturated by Am n \, (the con- 
tribution due to the explicit density dependence of the 
neutron-proton mass difference). Rather more important 

is Am n p , the contribution due to the difference in the 
wave functions of valence proton and neutron weighted 
by the local (density and density-gradient induced) vari- 
ation of the effective mass of the nucleon. In fact, when 
we restrict the in-medium reduction of the (averaged ef- 
fective) proton mass to about 1 % of the free proton mass, 
we obtain a rather precise description of the NSA. 



Of course, the calculations of AE^m in ref. [T7] (see 
eq. (fSTj) ) takes into account the effect due to the differ- 
ent wave functions of valence nucleonj 1 " 1 ! But in ref. [T7] 
a constant value of the nucleon mass was used in the cal- 
culations, namely the free mass. So it might be possible 
to resolve the anomaly of the mirror nuclei by consider- 
ing a dynamical (local) mass of the nucleon that is only 
slightly reduced in comparison to the free mass. In this 
respect, it is interesting to note that in the local-density- 
functional approach to many-body nuclear systems an ad- 
ditional term has been introduced into the Coulomb part 
of the energy density of the nuclear system that removes 
the anomaly [49,50, 51,-52] . Moreover, this term is chosen 
such that it is proportional to the isoscalar rather than the 
isovector densitM 12 !. In this connection, it is even argued 
in refs. [53|f54] that the effective-to- free-nucleon-mass ratio 
"t-n/ to n is unity to within a few percent. It seems that dif- 
ferent model calculations lead to similar conclusions about 
the origin of the Nolen-Schiffer anomaly. 

Returning to our model, we note that the results might 
be improved by a more self-consistent calculation that 
takes into account the feedback(s) between the modified 
skyrmion and the local nuclear background and by the 
inclusion of further degrees of freedom, which may any- 
how be needed from more detailed considerations about 
the nucleon structure and the nucleon-nucleon interaction. 
Also the non-local character of the effective nucleon mass 
may be of importance. Note that many-body calculations 
(sec e.g. 3,4,8,9,55]) indicate that, while m* is small in- 
side the Fermi sphere, it can reach values close to the 
vacuum value near k-p. In this connection, however, we 
should remark that the gradient terms which are present 
in our model do not noticeably affect the scaling behav- 
ior of m* . They are important for the surface behavior of 
Am* , though. 

We would like to thank Horst Lenske for useful discussions 
and for a computer code calculating nuclear densities of finite 
nuclei which allowed us to extract the background-density pa- 
rameters. We are also grateful to Frank Griimmer for provid- 

11 This phenomenon is known as Thomas-Ehrman effect [471 
EH]. 

12 Note, however, that the isoscalar contribution is surface 
dominated. 
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A Mass and moments of inertia of a classical 
soliton in the nonspherical scenario 

In the following, we use the definitions 

P r =d r P, Pg = d e P, r =d r 0, e = d e O, 

Sp = sin P, Cp = cos P, 5@ s sinS, 

C e = cos6>, S e = sind, C g = cosd. (A.l) 

The mass of the nonperturbated system of the Lagrange 
functional (l36l) has the form 



M£ P = it J drr 2 JsedO x {-^ 





)2 i 2 n2 



si 



p e z + r 2 p 7 z + s% ( -§ + e 2 e + r 2 el 



Sp 



§(PZ + r 2 P?) 



+ a^l(i + m -^(x))(i~cp) 



(A.2) 



where x = \lr 2 + 2rR cos 8 + R 2 [34] . For the moment of 
inertia, we introduce the following definition 



A = 2n J drr 2 J S g d8 A, 
o o 



(A.3) 



where the contributions of the different parts of the La- 
grange functional (|36[) are given as 



A* 



,12 ■ 



[l + C%) (P^ + r 2 P 2 ) 



2e 2 r 2 



s 2 P (^ + c 2 (0 2 + r 2 2 r \ 



(A.4) 



2e 2 r 2 



CeC ^(Pt + r 2 P? + Sl {&$+r 2 P?] 



r 2 P r (P r Og - O r Pg) + S%^0 



5/7 



(A.5) 



x * _F 2 f X f(x) 



S 2 P IS? 



P 2 +Sp { C 2 -M-+0 2 

•■II 



ff((l + C e 2 ) (Pt + S 2 P 2 ) 



2e 2 r 2 

2 / d2 i q2 /q2A s~i2 



{P 2 + S 2 P 2 r )C 2 )+r 2 {P r 0g-Pg0 r Y , (A.6) 



F: 



m 

P 2 + r 2 P 2 + S 2 P {0 2 + r 2 2 )) S 2 , (A.7) 



\* 2 77 I 1 i V*V \ q2 q2 

A Luf2,33~- ~ I 1 " 1 ) d0dp 

Ames " 
A* 



F 2 

— D Dp , 

F -fA*S 2 e S 2 P . 



(A.8) 
(A.9) 



B Charges and magnetic moments in the 
nonspherical scenario 

The zero component of the baryonic current in nonspher- 
ically deformed states has the form 

Bo = -r% ^| (P r 0o Po0 r ) ■ (B.10) 

In calculating the third component of the isospin current 
, one finds 

J d 3 rV/ (3) = (a*+u; 2 -n 3 )A^ 33 + A* cnv ^T 3 .(BAl) 

Consequently, the angle-averaged isospin charge density is 
given as 



T 3 = (a* +uj 2 - /2 3 )A* n , 33 + A*„ 
T 3 - A* c 



A 



A* 



uf2,33 ' A env i 



a;: 



(B.12) 



uif2,33 



where 



A cnv = 2 / Sin 6 d0X c^ 



sin0d#A^ 33 , (B.13) 
(B.14) 



Since the charges of the nucleons are defined as 



p%(r,6)d d r± P Ur,0)d 6 r 



(B.15) 
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the isoscalar and the isovector density distributions have 
here the following form 



18. B. A. Li, C. M. Ko and W. Bauer, Int. J. Mod. Phys. E 



n s - ^ - 
PE ~ 2 A 



A 



Kn,33 + Knv i (B.16) 



X 



uf2,33 



2A 



(B.17) 



Analogous calculations for the magnetic moments lead to 



*(5> 



P Ur,0)d 3 r± P l l (r,9)d 3 r, (B.18) 



where 



s _ m N (l + vi; nv ) 2 2 
Pm — 7TZ B r sm U i 



4A 



(B.19) 



a; 47, 33 



V _ TO N 



-f(l- x l{x))SlS% 



<? 2 



x (P 2 + r 2 P 2 + S£(0g + r 2 6> 2 )) S|l . (B.20) 
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